Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1578  769 (S20W36 L=305)   461    130        2 (1523/ 717-1618/ 804)  35.5 W31-W38 S23-S18
      2  192  651 (S18E72 L=197)  1428    113        3 ( 165/ 604- 219/ 710)  56.3 E78-E66 S20-S14
      3 1705  641 (S29W50 L=319)   448     54        2 (1677/ 616-1730/ 663)  -9.5 W46-W52 S30-S27
      4 1158  531 (S34W06 L=275)   560    123        1 (1102/ 503-1211/ 575)  14.0 W01-W09 S35-S30
      5 1383  497 (S37W23 L=292)  1035    279        0 (1281/ 441-1509/ 567)  22.4 W15-W32 S41-S32
      6 1015  508 (S34E05 L=264)   265     62        0 ( 982/ 490-1048/ 528) -36.9 E06-E02 S35-S32
      7 1093  310 (S50E01 L=268)   329     83        1 (1061/ 289-1137/ 326)   9.5 E03-W03 S51-S48
      8  864  201 (S58E29 L=240)   467     48        1 ( 837/ 181- 889/ 224) -21.8 E32-E26 S60-S55
  

    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