Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1670 1544 (N41W50 L=127)   935    173        2 (1611/1468-1731/1593) -48.8 W45-W54 N35-N44
      2 1072 1492 (N32E01 L= 76)   270     35        1 (1047/1472-1097/1512)  18.4 E02-W01 N30-N33
      3 1207 1412 (N28W09 L= 86)   390     62        1 (1186/1375-1241/1446) -90.0 W07-W11 N25-N29
      4  281 1291 (N12E57 L= 20)  1796    340        7 ( 249/1133- 316/1403) -80.0 E61-E52 N02-N19
      5  144 1246 (N07E76 L=  1)  1106    163        1 ( 128/1214- 163/1284) -38.7 E81-E71 N05-N09
      6  970  838 (S10E02 L= 75)   383    159        0 ( 953/ 759- 998/ 918) -71.6 E03-W00 S14-S05
      7 1751  724 (S12W56 L=132)   579    123        2 (1718/ 694-1795/ 749)  23.2 W52-W60 S13-S10
      8 1324  652 (S19W23 L= 99)   259    105        0 (1304/ 611-1350/ 693) -59.0 W21-W25 S21-S16
      9 1369  464 (S32W31 L=107)   396     88        0 (1325/ 449-1415/ 477)  18.4 W27-W33 S33-S30
  

    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