Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1323 1587 (N41W10 L=193)   529     97        0 (1291/1553-1358/1618) -31.0 W07-W13 N38-N43
      2  314 1532 (N17E68 L=114)   559     43        1 ( 299/1508- 332/1564) -45.0 E71-E65 N15-N19
      3 1657 1330 (N30W40 L=222)  1187    304        3 (1562/1215-1732/1413) -35.5 W31-W47 N23-N34
      4  370 1371 (N08E52 L=131)   410    132        1 ( 335/1337- 415/1406)  45.0 E54-E49 N06-N11
      5 1816 1168 (N22W56 L=239)   795     98        2 (1791/1119-1838/1211)  63.4 W53-W59 N19-N25
      6  183 1017 (S16E64 L=119)  2270    378        4 ( 114/ 882- 242/1162) -40.1 E80-E55 S23-S07
      7  199 1068 (S12E62 L=120)   518     62        3 ( 170/1045- 222/1092)  26.6 E66-E58 S13-S10
      8  687  769 (S22E17 L=166)  1891    474        1 ( 514/ 654- 856/ 911)  47.6 E29-E07 S33-S10
      9 1567  586 (S17W46 L=229)   303      9        0 (1545/ 568-1589/ 604)  14.0 W44-W48 S18-S16
     10 1543  509 (S22W48 L=231)   703    163        2 (1464/ 472-1605/ 550)  41.2 W43-W51 S26-S18
     11 1151  207 (S54W40 L=223)  1431    193        2 (1082/ 186-1250/ 238)  33.7 W31-W48 S57-S48
  

    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