Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  966 1699 (N51E01 L= 54)   415     82        0 ( 945/1661- 998/1728)  51.3 E03-W02 N47-N53
      2 1145 1572 (N40W13 L= 68)   559    159        1 (1090/1516-1208/1621) -48.8 W08-W16 N35-N43
      3  303 1443 (N32E62 L=354)  1170    151        1 ( 286/1383- 333/1509) -36.0 E68-E54 N28-N36
      4 1402 1486 (N32W31 L= 87)   469     96        1 (1350/1472-1456/1501) -12.5 W26-W35 N30-N32
      5 1831 1351 (N19W69 L=125)  1489    196        3 (1809/1261-1854/1437)  48.0 W64-W75 N13-N23
      6  997 1383 (N27E00 L= 55)   286     84        0 ( 959/1364-1041/1398)  14.0 E02-W02 N26-N28
      7 1699  787 (S15W47 L=102)  2215    345        3 (1597/ 644-1770/ 912)  50.9 W39-W52 S23-S07
      8  995  770 (S11E03 L= 53)   215     83        0 ( 955/ 743-1025/ 792)  31.0 E05-W00 S12-S09
      9 1327  684 (S18W18 L= 74)   252     60        0 (1305/ 660-1349/ 711)  71.6 W16-W19 S20-S16
     10 1093  542 (S27W02 L= 58)   484    174        0 (1028/ 515-1179/ 569)  16.7 E02-W08 S28-S25
     11 1182  164 (S65W14 L= 70)   751     81        0 (1149/ 143-1210/ 187) -26.6 W08-W20 S68-S60
  

    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