Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1118 1835 (N48E24 L=120)   436    107        1 (1073/1808-1167/1852)  71.6 E25-E21 N43-N52
      2  571 1606 (N16E46 L= 98)   690    122        2 ( 552/1539- 596/1658) -41.2 E50-E41 N11-N19
      3  695 1587 (N17E36 L=107)   607    190        2 ( 624/1565- 783/1625)  29.1 E41-E30 N14-N19
      4 1595 1520 (N40W25 L=168)   563    131        0 (1545/1467-1630/1571) -26.6 W20-W28 N37-N42
      5  749 1480 (N12E28 L=115)   265    125        0 ( 736/1421- 766/1536) -54.5 E30-E25 N08-N15
      6  710  888 (S23E15 L=128)   260    127        0 ( 683/ 831- 730/ 951)  90.0 E15-E14 S26-S18
      7  546  727 (S36E24 L=120)   566    186        1 ( 468/ 693- 624/ 769)  33.7 E30-E18 S39-S31
      8  983  738 (S24W05 L=149)   186     48        0 ( 962/ 725-1014/ 758)   0.0 W03-W07 S25-S23
      9 1227  687 (S20W21 L=165)   298    110        0 (1206/ 628-1242/ 735) -39.8 W18-W24 S23-S17
     10  751  639 (S36E06 L=137)   380    159        1 ( 730/ 575- 793/ 709) -45.0 E09-E01 S39-S32
     11 1539  306 (S30W69 L=213)   896    110        4 (1483/ 288-1599/ 326)  13.2 W61-W78 S31-S27
     12  832  223 (S65W24 L=167)  1480    133        2 ( 784/ 176- 882/ 281) -19.3 W14-W35 S68-S60
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y