Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1288 1634 (N32W38 L=100)   599    143        1 (1234/1595-1350/1673)   0.0 W32-W44 N30-N33
      2 1377 1464 (N20W35 L= 97)   288     63        0 (1360/1434-1395/1501)  36.9 W32-W37 N18-N21
      3 1852 1328 (S03W68 L=130)   587     91        1 (1829/1287-1876/1369)  68.2 W66-W70 S06-N00
      4 1242 1268 (N13W19 L= 81)   295    103        0 (1188/1247-1290/1296) -18.4 W16-W22 N12-N15
      5  751 1160 (N21E12 L= 50)  1823    472        3 ( 565/1088- 987/1235) -36.5 E23-W00 N11-N29
      6 1721  964 (S18W41 L=104)  1721    420        3 (1685/ 773-1743/1164)  55.0 W34-W48 S28-S07
      7 1808 1107 (S13W53 L=115)   656     58        1 (1780/1078-1833/1135) -36.9 W50-W55 S14-S10
      8 1298 1092 (N02W17 L= 79)   318     52        1 (1275/1065-1324/1118)  18.4 W15-W18 N01-N04
      9  819  737 (S04E19 L= 43)   267     67        0 ( 790/ 696- 856/ 768) -80.5 E20-E17 S07-S01
     10 1481  646 (S29W17 L= 80)  1441    276        1 (1372/ 570-1546/ 762)  24.8 W09-W23 S32-S23
     11  957  691 (S11E13 L= 50)   448    147        0 ( 875/ 659-1032/ 734) -11.3 E17-E07 S11-S08
     12 1038  603 (S18E11 L= 51)   288    154        0 (1014/ 533-1063/ 671)  45.0 E14-E07 S21-S14
  

    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