Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  612 1817 (N41E55 L= 93)  1340     94        3 ( 586/1778- 639/1859) -19.3 E67-E47 N37-N45
      2  963 1673 (N34E16 L=132)   403     76        1 ( 930/1649-1000/1694)  36.9 E17-E13 N32-N36
      3  384 1605 (N23E59 L= 88)  1267    165        4 ( 356/1523- 415/1676) -32.3 E71-E51 N16-N28
      4  917 1384 (N14E12 L=136)   417     85        3 ( 879/1355- 957/1415)  59.0 E13-E09 N11-N16
      5 1808 1261 (N22W53 L=200)  2735    300        3 (1732/1145-1880/1380) -36.3 W45-W61 N15-N29
      6 1519 1312 (N19W27 L=175)   312     39        0 (1499/1283-1537/1343)  78.7 W25-W28 N16-N21
      7 1105 1248 (N08W02 L=149)   355    131        1 (1059/1218-1169/1282)  -9.5 E01-W05 N06-N09
      8  308 1241 (S01E51 L= 96)   332     62        0 ( 295/1207- 321/1276)  90.0 E52-E49 S03-N00
      9  150 1125 (S09E68 L= 80)   815    121        3 ( 129/1054- 175/1174) -53.1 E70-E63 S13-S05
     10 1268 1095 (N02W13 L=161)   347    135        1 (1223/1048-1306/1120)  53.1 W11-W15 S01-N03
     11 1369 1064 (N01W20 L=168)   288    110        0 (1322/1039-1406/1106)  45.0 W17-W21 N00-N04
     12 1312  962 (S05W18 L=166)   555    236        1 (1211/ 909-1420/ 992)  -4.1 W11-W25 S06-S04
     13 1344  783 (S15W23 L=171)   802    254        0 (1284/ 682-1411/ 874)  77.0 W21-W25 S21-S08
     14 1588  699 (S16W43 L=190)   274     71        0 (1567/ 680-1625/ 728)  -9.5 W40-W46 S16-S14
     15 1388  636 (S23W30 L=178)   386     94        1 (1342/ 616-1430/ 662)  -8.1 W26-W33 S24-S21
     16 1808  578 (S17W71 L=219)   978     82        0 (1766/ 557-1851/ 596)  17.1 W66-W80 S19-S15
     17  442  485 (S48E44 L=104)   465     80        1 ( 410/ 447- 470/ 510) -59.0 E46-E40 S51-S46
     18 1608  444 (S30W58 L=206)   414     52        0 (1579/ 424-1629/ 465)  33.7 W54-W60 S30-S28
  

    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