Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1387 1776 (N37W56 L=166)   786    123        2 (1317/1749-1445/1806) -13.0 W49-W63 N34-N38
      2 1066 1614 (N37W22 L=132)   360     58        1 (1037/1585-1087/1644)  31.0 W19-W24 N35-N38
      3 1045 1380 (N24W11 L=121)   317     75        1 (1008/1357-1087/1397) -36.9 W09-W13 N22-N25
      4 1839 1337 (S01W65 L=175)  1276    106        3 (1791/1301-1892/1370) -18.4 W59-W73 S03-N01
      5 1288 1153 (N05W18 L=128)   196     96        0 (1242/1128-1327/1179)   0.0 W14-W20 N04-N05
      6  501 1130 (N23E29 L= 81)  1477    320        4 ( 336/1092- 639/1177) -24.4 E42-E20 N19-N29
      7 1068 1129 (N10W05 L=115)   191    110        0 (1017/1118-1123/1141) -18.4 W02-W08 N08-N11
      8 1504  867 (S16W24 L=134)  1467    437        3 (1348/ 731-1613/ 994)  22.8 W11-W31 S19-S09
      9  781  902 (N04E16 L= 94)   460    177        0 ( 715/ 878- 836/ 946)  -8.1 E19-E12 N02-N07
     10  153  856 (N14E68 L= 42)  1241    119        3 (  93/ 832- 214/ 890) -11.9 E79-E60 N12-N16
     11  353  869 (N12E46 L= 64)   418     85        1 ( 322/ 843- 391/ 890)   0.0 E49-E42 N10-N13
     12 1146  878 (S06W03 L=113)   191     86        0 (1103/ 857-1190/ 890) -31.0 W00-W05 S08-S04
     13  372  819 (N08E45 L= 65)   253      0        0 ( 354/ 783- 394/ 852)  36.9 E47-E43 N06-N10
     14 1724  704 (S33W39 L=149)  1760    305        3 (1652/ 608-1817/ 844)  30.5 W31-W49 S38-S27
     15  318  644 (S00E56 L= 54)  1821    324        2 ( 227/ 526- 409/ 791)  20.0 E77-E43 S06-N06
     16 1491  677 (S27W19 L=129)   261     77        0 (1466/ 638-1514/ 714)  31.0 W16-W21 S28-S25
     17  632  636 (S07E33 L= 77)   431    123        0 ( 592/ 589- 682/ 677)  14.0 E36-E28 S09-S05
     18  951  220 (S43E35 L= 76)  1392    354        3 ( 806/ 169-1083/ 318)   5.4 E50-E18 S46-S39
  

    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