Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  445 1646 (N23E62 L= 41)   508     81        1 ( 408/1624- 483/1661)   0.0 E67-E55 N22-N24
      2 1122 1525 (N30E06 L= 97)   265     62        0 (1087/1504-1147/1544) -11.3 E08-E03 N29-N31
      3  869 1404 (N17E18 L= 85)   453    123        1 ( 815/1365- 919/1431)  15.9 E21-E14 N14-N18
      4 1149 1370 (N21E00 L=103)   427    116        0 (1112/1342-1181/1406)  33.7 E02-W01 N19-N22
      5 1157 1276 (N15W02 L=105)   705    161        4 (1099/1196-1195/1328) -39.8 E01-W05 N10-N18
      6  201 1061 (S16E61 L= 41)   558     61        0 ( 180/1029- 218/1101)  78.7 E64-E59 S18-S13
      7 1789 1008 (N14W53 L=156)   447    100        0 (1776/ 954-1808/1056) -51.3 W50-W55 N11-N16
      8  813  655 (S29E05 L= 98)  1142    292        0 ( 674/ 594- 938/ 721)  42.9 E14-W01 S35-S21
      9 1638  423 (S23W66 L=168)  1131     70        0 (1599/ 394-1667/ 456)  51.3 W61-W70 S25-S20
     10  710  396 (S50E08 L= 95)   493     89        1 ( 677/ 373- 761/ 418)  40.6 E10-E03 S53-S47
     11  988  332 (S47W19 L=121)  1476    263        0 ( 871/ 313-1110/ 351)  25.2 W08-W27 S51-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