Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1238 1454 (N27W24 L=183)   546    111        2 (1194/1415-1289/1487)  12.5 W19-W28 N25-N29
      2  450 1254 (N31E35 L=124)  1576    441        1 ( 349/1104- 583/1396) -49.2 E45-E26 N20-N42
      3  188 1213 (N32E60 L= 99)   479     70        1 ( 167/1191- 215/1231) -11.3 E62-E57 N30-N32
      4  197 1117 (N26E59 L=100)   924    157        3 ( 170/1078- 231/1154)  53.1 E62-E55 N23-N28
      5 1648 1088 (S03W41 L=200)   786    158        3 (1586/1062-1700/1118) -18.4 W35-W44 S05-S01
      6  637 1063 (N16E23 L=136)   329    118        1 ( 593/1047- 699/1082) -15.9 E26-E19 N15-N18
      7 1490 1019 (S03W28 L=187)   353    199        3 (1471/ 964-1507/1076)  56.3 W26-W30 S06-N00
      8  260  992 (N18E54 L=105)   371    113        0 ( 244/ 943- 270/1059)  66.8 E56-E52 N14-N21
      9 1630  976 (S09W37 L=196)   446    220        2 (1566/ 901-1697/1024)  14.0 W30-W42 S12-S07
     10  804  906 (N04E15 L=144)   998    256        2 ( 723/ 825- 884/1012) -63.4 E18-E11 S01-N11
     11 1513  818 (S16W26 L=185)   427     59        0 (1498/ 779-1529/ 861)  78.7 W24-W27 S18-S13
     12  799  570 (S16E24 L=136)   866    188        3 ( 700/ 504- 853/ 607) -48.8 E28-E19 S20-S11
     13  662  276 (S35E47 L=112)   929    183        1 ( 625/ 234- 690/ 316)  28.6 E53-E42 S37-S31
  

    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