Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  995 1549 (N24E17 L= 15)   234     58        0 ( 967/1526-1029/1566)  53.1 E18-E15 N21-N25
      2 1112 1500 (N24E09 L= 24)   459    116        0 (1061/1461-1166/1528)   0.0 E12-E04 N22-N25
      3 1669 1463 (N41W33 L= 65)   409     51        0 (1641/1444-1692/1488)  18.4 W29-W35 N39-N43
      4 1872 1122 (N27W58 L= 91)   658    102        1 (1852/1068-1890/1170) -51.3 W55-W61 N23-N29
      5 1652 1082 (N16W36 L= 69)   615    211        2 (1604/1005-1677/1157) -49.4 W31-W38 N11-N19
      6  623 1011 (S17E23 L= 10)   388    124        1 ( 603/ 947- 655/1057)  90.0 E24-E21 S20-S13
      7 1783  960 (N14W51 L= 83)   374    106        1 (1757/ 934-1823/ 996)  -7.1 W47-W55 N11-N14
      8  744  919 (S19E13 L= 19)   253     82        0 ( 698/ 905- 787/ 937)  31.0 E15-E10 S20-S17
      9  304  781 (S39E47 L=346)   499     57        1 ( 263/ 764- 341/ 795)   5.7 E52-E42 S39-S37
     10 1560  626 (S11W45 L= 78)   308     43        1 (1531/ 608-1583/ 650)  26.6 W43-W47 S12-S10
  

    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