Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1351 1570 (N35W04 L=114)  2041    404        2 (1212/1378-1447/1710) -33.1 E09-W15 N24-N41
      2  282 1526 (N06E69 L= 40)  2013    299        4 ( 247/1444- 332/1599)  90.0 E76-E65 N00-N11
      3  613 1559 (N13E40 L= 69)   633    146        1 ( 546/1530- 686/1599)  40.6 E44-E36 N10-N17
      4  930 1468 (N16E17 L= 92)   373    117        0 ( 874/1435- 986/1487)   7.1 E21-E13 N14-N16
      5 1016 1283 (N07E07 L=102)   289     86        0 ( 988/1246-1047/1324) -31.0 E09-E04 N05-N08
      6  139 1258 (S12E73 L= 36)   740     70        0 ( 125/1216- 156/1287)  90.0 E75-E71 S14-S10
      7  265 1120 (S19E53 L= 57)   579    127        1 ( 248/1053- 290/1183)  90.0 E54-E51 S22-S14
      8 1337  595 (S21W31 L=141)   621    160        1 (1256/ 579-1416/ 617)  29.1 W26-W35 S23-S18
      9 1501  385 (S27W54 L=163)   410     59        1 (1473/ 370-1534/ 399)  21.8 W51-W56 S28-S25
  

    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