Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  437 1428 (N33E36 L=190)   560     92        1 ( 396/1394- 474/1463)  31.0 E39-E32 N31-N35
      2  615 1251 (N19E23 L=203)   349    101        0 ( 585/1193- 650/1299)  50.2 E24-E19 N15-N21
      3 1169 1269 (N12W12 L=238)   193     53        0 (1144/1240-1196/1290) -45.0 W10-W13 N10-N13
      4  530  888 (S01E32 L=194)   271     90        0 ( 487/ 851- 561/ 913) -36.9 E34-E30 S03-N00
      5  279  766 (S05E55 L=171)   442    113        0 ( 250/ 716- 297/ 811) -90.0 E56-E52 S08-S02
      6  282  656 (S12E58 L=168)   876    311        1 ( 238/ 579- 325/ 744)  31.0 E66-E51 S15-S06
      7  464  565 (S20E45 L=181)   807    265        3 ( 414/ 490- 522/ 654)  28.1 E52-E37 S23-S14
      8 1415  324 (S53W23 L=248)   318     55        2 (1391/ 311-1444/ 342) -23.2 W19-W26 S55-S51
  

    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