Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1550 1700 (N32W60 L= 25)   731    144        1 (1519/1671-1586/1718)   0.0 W55-W63 N30-N33
      2 1015 1584 (N41W15 L=340)   448     76        1 ( 986/1545-1036/1620)  59.0 W12-W17 N38-N43
      3 1752 1416 (N11W62 L= 27)   700    125        2 (1723/1368-1794/1461) -71.6 W58-W64 N07-N13
      4 1004 1421 (N31W08 L=334)   397    111        0 ( 956/1400-1062/1448) -39.8 W05-W11 N28-N33
      5 1455 1409 (N19W37 L=  2)   286     78        0 (1419/1394-1491/1422) -26.6 W34-W39 N17-N20
      6 1635 1388 (N13W49 L= 15)   322     83        0 (1596/1373-1669/1399) -18.4 W46-W52 N11-N14
      7 1307 1140 (N08W19 L=345)   885    268        2 (1261/1003-1351/1268)  82.9 W16-W22 S01-N15
      8  954 1104 (N13E02 L=323)   216     88        0 ( 926/1061- 982/1152) -71.6 E03-E01 N10-N16
      9  271  848 (N09E56 L=270)   426    106        0 ( 251/ 820- 305/ 900) -68.2 E57-E52 N07-N12
     10 1427  716 (S20W18 L=344)   454    129        0 (1384/ 653-1474/ 769)  35.5 W14-W22 S23-S17
     11 1216  408 (S36E02 L=324)   405     88        2 (1189/ 379-1245/ 447)  14.0 E04-W00 S38-S32
  

    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