Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  708 1729 (N31E43 L=352)  1127    187        0 ( 664/1672- 766/1785) -19.7 E50-E36 N28-N34
      2  660 1555 (N19E37 L=358)  2090    355        3 ( 516/1388- 726/1634) -28.8 E49-E29 N08-N23
      3 1206 1478 (N28E01 L= 35)  1367    446        3 (1031/1338-1402/1611) -13.1 E15-W15 N23-N33
      4  484 1455 (N09E46 L=350)   828    221        3 ( 447/1381- 523/1526) -53.1 E49-E40 N05-N13
      5 1618 1429 (N37W30 L= 66)   539    133        0 (1560/1382-1672/1466) -11.3 W25-W35 N35-N38
      6  832 1363 (N11E20 L= 16)   324     87        0 ( 803/1319- 865/1405) -33.7 E22-E16 N09-N13
      7 1758 1288 (N32W44 L= 79)   533    103        0 (1716/1247-1786/1312)  -9.5 W39-W47 N30-N33
      8  891 1167 (N01E11 L= 24)   245     88        1 ( 852/1139- 932/1191)  53.1 E13-E09 S01-N03
      9 1738 1051 (N16W45 L= 81)   310    102        1 (1710/1017-1757/1100)  78.7 W43-W47 N14-N19
     10  884  491 (S39W06 L= 42)   380     78        0 ( 850/ 465- 931/ 523)   8.1 W03-W10 S39-S37
  

    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