Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1883 (N56E42 L= 36)  1607    213        2 ( 864/1851-1042/1910)  54.0 E48-E36 N50-N62
      2  665 1570 (N23E37 L= 40)   416     90        1 ( 624/1548- 713/1591)   0.0 E41-E33 N22-N24
      3 1286 1271 (N18W11 L= 88)  1172    345        6 (1119/1205-1420/1327)   3.0 W00-W19 N15-N20
      4 1493 1249 (N21W25 L=102)   372     85        1 (1467/1212-1514/1291)  90.0 W23-W26 N18-N24
      5 1508 1172 (N17W27 L=104)   366     94        2 (1470/1143-1542/1194)  31.0 W24-W29 N15-N18
      6 1350 1154 (N13W17 L= 94)   198     73        0 (1325/1121-1371/1185) -26.6 W14-W18 N10-N13
      7 1605  603 (S14W49 L=127)   386     42        2 (1578/ 585-1634/ 622)   0.0 W46-W52 S14-S12
      8 1203  554 (S26W23 L=100)   864    205        1 (1119/ 524-1306/ 611)  38.7 W17-W27 S29-S20
      9  975  488 (S36W10 L= 87)   360    109        0 ( 928/ 476-1030/ 506)  15.9 W06-W13 S37-S33
     10  948  315 (S50W16 L= 93)   333     54        0 ( 926/ 291- 969/ 340) -18.4 W12-W18 S51-S48
     11 1161  256 (S48W38 L=115)   743     81        1 (1125/ 231-1191/ 301) -60.3 W34-W42 S50-S43
     12 1268  242 (S46W49 L=126)  1596    165        1 (1200/ 216-1347/ 272)   6.3 W39-W57 S48-S42
  

    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