Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1392 1678 (N45W04 L= 14)  2269    370        2 (1276/1551-1525/1791) -23.7 E08-W17 N39-N52
      2  470 1722 (N21E66 L=305)  1291    134        5 ( 410/1671- 512/1765)  90.0 E68-E62 N16-N24
      3  725 1663 (N22E39 L=331)  1249    256        3 ( 607/1594- 806/1721)  63.4 E43-E35 N15-N28
      4 1201 1584 (N31E06 L=  4)   464    110        1 (1173/1526-1219/1633) -54.5 E09-E04 N26-N33
      5  861 1551 (N19E25 L=345)   268    108        0 ( 832/1535- 909/1587)  -9.5 E27-E21 N17-N20
      6 1434 1442 (N29W13 L= 23)   442    197        0 (1356/1380-1500/1479)   0.0 W07-W18 N27-N30
      7  874 1260 (N02E15 L=356)   249     46        0 ( 837/1240- 899/1275)  14.0 E16-E12 N01-N03
      8 1876 1169 (N29W57 L= 67)   614     77        1 (1851/1134-1903/1208) -23.2 W53-W62 N26-N30
      9  199 1127 (S20E62 L=308)  1102    143        6 ( 173/1058- 221/1197)  83.7 E64-E59 S24-S15
     10  313  993 (S26E47 L=323)   405    106        0 ( 269/ 959- 351/1035)  40.6 E51-E43 S28-S22
     11  839  851 (S21E06 L=  4)   300    113        0 ( 781/ 840- 893/ 863)  23.2 E09-E02 S23-S20
     12  241  787 (S40E56 L=314)   692    118        1 ( 195/ 743- 294/ 816)  13.2 E65-E48 S43-S38
     13 1431  388 (S29W51 L= 61)  1318    229        1 (1369/ 301-1520/ 448) -81.3 W47-W56 S35-S22
  

    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