Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1096 1507 (N24E10 L= 24)   484    121        1 (1046/1468-1153/1541)   0.0 E13-E05 N22-N25
      2 1635 1485 (N42W29 L= 63)   784    131        4 (1581/1439-1680/1529) -11.3 W23-W34 N38-N44
      3 1636 1101 (N17W35 L= 68)   563    230        2 (1577/1011-1666/1168) -58.0 W29-W37 N11-N19
      4 1863 1117 (N26W57 L= 91)   551     81        1 (1843/1074-1883/1157) -26.6 W54-W60 N24-N28
      5  612 1025 (S16E24 L=  9)   348    110        1 ( 592/ 964- 641/1069)  90.0 E25-E23 S19-S13
      6 1770  963 (N13W49 L= 83)   332     79        1 (1751/ 937-1803/1000)   0.0 W46-W52 N11-N14
      7  730  927 (S19E14 L= 19)   246    118        0 ( 677/ 915- 785/ 940)  33.7 E17-E11 S20-S16
      8  299  786 (S39E48 L=346)   484     77        2 ( 260/ 771- 335/ 811)  12.5 E52-E43 S39-S36
      9 1544  633 (S11W44 L= 77)   597    123        2 (1496/ 609-1594/ 657)   6.3 W39-W48 S13-S09
  

    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