Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  732 1763 (N30E47 L=309)   493    121        2 ( 709/1704- 756/1808) -45.0 E51-E43 N25-N33
      2 1220 1698 (N41E11 L=345)   607    149        1 (1180/1661-1256/1754) -26.6 E16-E06 N38-N44
      3  769 1510 (N14E29 L=326)   742    260        1 ( 738/1386- 791/1627) -51.8 E35-E23 N07-N22
      4  884 1540 (N19E24 L=332)   708    177        5 ( 846/1448- 910/1626) -42.0 E29-E19 N14-N24
      5  570 1465 (N07E40 L=316)   328    143        0 ( 545/1408- 607/1522) -82.9 E41-E39 N03-N11
      6  317 1452 (N01E61 L=295)   663    112        1 ( 268/1429- 371/1474)  18.4 E66-E56 S01-N02
      7 1688 1374 (N36W34 L= 30)  1156    162        4 (1655/1294-1721/1448) -56.3 W30-W38 N30-N40
      8 1003 1300 (N09E09 L=347)  1057    244        3 ( 922/1226-1103/1375)  65.2 E11-E05 N02-N15
      9  360 1193 (S13E47 L=309)   521    153        0 ( 299/1153- 425/1225)  -4.8 E52-E40 S14-S12
     10 1220 1060 (N01W10 L=  5)   962    262        2 (1139/1023-1302/1134)   0.0 W03-W14 S01-N05
     11  763  974 (S16E13 L=342)   779    274        3 ( 637/ 917- 870/1032)  34.7 E21-E08 S20-S09
     12 1678  925 (N08W42 L= 38)   505    102        2 (1663/ 867-1693/ 980) -45.0 W39-W45 N05-N11
     13  278  802 (S39E50 L=306)  1179    206        2 ( 195/ 761- 376/ 858)  20.6 E63-E39 S42-S33
     14 1100  607 (S27W18 L= 13)  1743    533        4 ( 891/ 531-1336/ 698)  45.0 W06-W27 S37-S14
     15  414  527 (S55E36 L=320)   946    120        2 ( 365/ 500- 463/ 564)   4.1 E43-E29 S57-S51
     16 1516  416 (S23W56 L= 51)   374     82        2 (1474/ 403-1560/ 429)  33.7 W52-W59 S25-S21
     17 1397  274 (S35W61 L= 56)  1727    173        1 (1333/ 212-1467/ 331) -70.0 W57-W66 S40-S29
     18 1161  248 (S46W46 L= 42)   455     59        2 (1144/ 218-1182/ 287) -18.4 W41-W50 S48-S44
  

    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