Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1027 1760 (N38E24 L=100)  1095    185        2 ( 968/1701-1122/1792)  71.6 E27-E19 N31-N43
      2 1487 1497 (N34W16 L=140)   461     51        0 (1466/1463-1512/1530) -36.9 W13-W17 N31-N36
      3  939 1333 (N08E13 L=112)   669    185        3 ( 845/1303-1030/1357)  26.6 E18-E07 N05-N10
      4 1671 1304 (N28W33 L=158)   518    106        0 (1637/1255-1699/1349) -29.7 W29-W36 N25-N29
      5  471 1217 (S09E38 L= 87)   368    129        0 ( 431/1158- 500/1275)  76.0 E39-E36 S13-S05
      6  368 1179 (S13E45 L= 80)   403     74        3 ( 346/1136- 392/1212)  68.2 E47-E42 S16-S11
      7  881 1090 (S07E09 L=115)   764    160        2 ( 823/1037- 927/1163) -35.5 E13-E05 S10-S03
      8  625 1110 (S12E25 L=100)   347     70        1 ( 599/1082- 653/1145)  63.4 E26-E23 S13-S09
      9 1888 1064 (N21W59 L=184)   503     98        1 (1865/1019-1911/1109) -29.7 W55-W63 N19-N23
     10 1808  956 (N12W52 L=177)   608    152        0 (1760/ 923-1856/ 982)  32.0 W48-W57 N09-N14
     11  930  912 (S15E02 L=122)   547    230        3 ( 877/ 872-1025/ 949)  49.4 E06-W02 S18-S11
     12 1042  888 (S14W05 L=129)   211     96        0 ( 996/ 876-1087/ 900)  21.8 W02-W07 S15-S12
     13 1785  860 (N06W53 L=177)   508     77        1 (1764/ 829-1816/ 898) -76.0 W50-W55 N04-N08
     14 1825  782 (N04W60 L=184)   593     46        1 (1799/ 760-1855/ 812)  31.0 W57-W62 N02-N05
     15 1101  708 (S22W13 L=138)   198    100        0 (1056/ 695-1141/ 726)  45.0 W10-W15 S23-S19
     16 1638  598 (S12W52 L=176)   458    102        2 (1611/ 556-1677/ 650)  90.0 W49-W53 S14-S07
     17 1675  491 (S16W62 L=187)   817     70        3 (1635/ 466-1712/ 515)  23.2 W57-W67 S17-S13
     18 1413  385 (S30W49 L=174)   560     57        1 (1377/ 362-1444/ 419)  21.8 W46-W52 S32-S27
  

    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