Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  736 1733 (N32E42 L=351)  1608    315        3 ( 688/1661- 798/1826) -30.5 E50-E33 N28-N39
      2  695 1582 (N21E36 L=357)  1325    297        5 ( 581/1499- 761/1623)  -7.1 E44-E28 N17-N24
      3 1160 1535 (N30E05 L= 27)   605    185        0 (1075/1458-1245/1594) -15.9 E12-W02 N27-N32
      4  504 1442 (N08E43 L=349)   588    163        1 ( 475/1395- 560/1497) -24.0 E47-E38 N05-N11
      5  729 1430 (N13E28 L=  5)   479    125        2 ( 700/1362- 755/1488) -76.0 E29-E26 N08-N16
      6 1639 1421 (N37W32 L= 65)   624    134        1 (1585/1368-1699/1456) -10.3 W27-W38 N35-N38
      7  869 1353 (N11E17 L= 16)   352    113        0 ( 838/1297- 899/1400) -29.7 E20-E13 N09-N13
      8 1379 1353 (N25W13 L= 46)   384    114        0 (1324/1315-1442/1376)   6.3 W09-W18 N23-N25
      9 1779 1282 (N32W46 L= 79)   460     93        1 (1741/1254-1810/1300)   0.0 W42-W49 N31-N33
     10  925 1153 (N01E09 L= 24)   205    112        0 ( 859/1140- 973/1173)  26.6 E12-E06 N00-N03
     11  909  479 (S39W09 L= 41)   408     66        0 ( 876/ 452- 947/ 512)  -8.1 W05-W12 S40-S36
  

    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