Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1710 (N36E24 L=196)   608     78        3 ( 819/1683- 901/1736)  36.9 E25-E20 N33-N37
      2  953 1598 (N29E13 L=207)   275     72        0 ( 926/1561- 985/1618)  68.2 E14-E11 N25-N30
      3  754 1207 (N01E18 L=201)   726    243        0 ( 717/1089- 792/1327) -58.4 E22-E14 S05-N08
      4 1611 1214 (N14W34 L=254)   307     56        1 (1585/1189-1637/1237) -14.0 W32-W36 N12-N15
      5 1907 1150 (N16W63 L=282)  1896    233        1 (1855/1099-1944/1195) -45.0 W56-W68 N13-N19
      6  572  981 (S14E27 L=192)  1058    306        3 ( 461/ 867- 692/1090)  45.0 E34-E20 S21-S06
      7 1909 1056 (N11W64 L=284)   764    152        0 (1863/1035-1936/1087)  29.1 W58-W68 N08-N13
      8  376  991 (S15E42 L=178)   832    277        4 ( 275/ 931- 481/1038)  20.6 E50-E34 S18-S11
      9 1207  811 (S16W14 L=234)   467    198        1 (1140/ 757-1244/ 859) -26.6 W09-W17 S19-S13
     10 1259  734 (S20W19 L=238)   197     77        0 (1209/ 720-1297/ 750)   9.5 W15-W21 S21-S18
     11  494  653 (S35E33 L=187)   433     48        1 ( 466/ 629- 524/ 682) -45.0 E35-E30 S36-S32
     12  861  215 (S65W02 L=222)   510     68        1 ( 834/ 191- 889/ 244)  -8.1 E01-W06 S67-S61
  

    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