Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  406 1206 (N19E38 L=310)  1818    291        3 ( 336/1079- 476/1324) -78.7 E43-E32 N11-N26
      2  914 1271 (N16E03 L=345)   339     87        0 ( 867/1245- 963/1307) -38.7 E05-W00 N13-N18
      3  260 1053 (N11E51 L=297)   490    102        2 ( 240/1006- 280/1111)  71.6 E53-E49 N08-N14
      4  837  778 (S12E15 L=334)   248    103        0 ( 825/ 724- 847/ 833)  71.6 E15-E13 S15-S09
      5 1450  652 (S28W23 L= 11)  1046    218        1 (1400/ 572-1480/ 753) -90.0 W19-W25 S33-S21
      6  925  686 (S19E11 L=337)   280     81        0 ( 880/ 662- 972/ 723)  18.4 E13-E07 S20-S17
      7  408  569 (S19E49 L=299)  1405    309        0 ( 376/ 475- 444/ 693)  42.7 E56-E42 S24-S12
      8 1216  506 (S35W06 L=354)   405    121        2 (1191/ 456-1241/ 564) -90.0 W03-W07 S38-S30
  

    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