Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  278 1536 (N16E71 L=349)   772     60        0 ( 258/1510- 297/1564)  53.1 E73-E68 N14-N18
      2  548 1483 (N17E41 L= 19)   262    114        1 ( 521/1425- 580/1535) -82.9 E41-E40 N13-N21
      3  790 1454 (N20E24 L= 37)   382    118        1 ( 739/1436- 847/1477)  23.2 E26-E19 N18-N21
      4  745 1177 (N03E20 L= 41)   254    104        2 ( 727/1134- 776/1230) -59.0 E20-E16 N00-N05
      5 1131  814 (S12W10 L= 71)   345    120        1 (1099/ 766-1146/ 876) -63.4 W08-W11 S14-S08
      6  984  808 (S15W02 L= 62)   191     82        1 ( 938/ 784-1019/ 836)  45.0 W00-W03 S16-S12
      7  865  683 (S25E03 L= 57)   279    130        0 ( 795/ 662- 926/ 696)  15.9 E07-W00 S26-S23
      8 1723  520 (S17W64 L=124)  1691    207        2 (1650/ 436-1773/ 579)  65.6 W59-W68 S23-S12
      9 1329  485 (S28W33 L= 93)  2419    664        7 (1081/ 396-1613/ 547)  27.2 W19-W54 S39-S20
  

    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