Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1410 (N25E13 L= 50)   763    196        0 ( 665/1366- 835/1437) -22.6 E19-E07 N22-N28
      2  220 1250 (N22E58 L=  4)  1707    256        2 ( 149/1154- 307/1335) -27.5 E73-E47 N14-N27
      3  677 1208 (N14E19 L= 43)   196     61        2 ( 650/1183- 705/1236) -33.7 E21-E17 N12-N15
      4  637 1020 (N03E23 L= 39)   234     50        0 ( 610/1003- 667/1038) -26.6 E25-E21 N01-N04
      5  696  977 (S00E20 L= 42)   413    201        1 ( 626/ 952- 786/1023)   5.7 E24-E14 S01-N01
      6  934  852 (S10E07 L= 55)   278     57        0 ( 903/ 825- 966/ 875)   0.0 E09-E05 S11-S08
      7 1650  724 (S26W41 L=103)   919    280        0 (1624/ 617-1681/ 855) -82.4 W37-W45 S32-S17
      8 1739  656 (S31W51 L=113)   779     91        1 (1696/ 635-1786/ 676)   6.3 W46-W57 S32-S29
      9  362  497 (S25E60 L=  3)  1854    212        3 ( 334/ 385- 399/ 579)  29.7 E76-E52 S32-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