Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1773 (N58W13 L= 60)   467     70        1 ( 932/1749-1007/1791)   0.0 W09-W18 N56-N59
      2  862 1755 (N59W01 L= 48)   327     52        1 ( 841/1729- 886/1774)  45.0 E01-W04 N56-N60
      3  771 1734 (N59E10 L= 37)   304     38        1 ( 747/1715- 799/1748)  -8.1 E13-E06 N57-N60
      4 1689 1537 (N25W63 L=110)   421     69        1 (1661/1511-1718/1567) -45.0 W60-W65 N22-N26
      5 1194 1487 (N32W20 L= 66)   380    186        0 (1102/1434-1249/1511) -37.9 W13-W22 N27-N34
      6  247 1331 (N33E61 L=346)  3688    434        7 ( 172/1130- 337/1491)  71.6 E66-E53 N20-N44
      7  402 1430 (N39E45 L=  2)   332     91        0 ( 374/1385- 425/1473) -41.2 E49-E41 N35-N42
      8 1007 1413 (N30W05 L= 52)   367    157        0 ( 953/1381-1077/1460) -45.0 W02-W09 N27-N34
      9 1687 1330 (N13W51 L= 98)   510    104        0 (1670/1294-1714/1387)  80.5 W48-W53 N10-N16
     10  365 1202 (N25E44 L=  2)   623     89        1 ( 332/1152- 404/1245) -40.6 E48-E40 N21-N27
     11 1852  940 (S13W59 L=106)   743     84        3 (1828/ 904-1882/ 976) -51.3 W56-W63 S15-S10
     12 1395  783 (S14W19 L= 66)   428     91        2 (1354/ 752-1434/ 808)   9.5 W16-W22 S15-S11
     13  525  579 (S15E40 L=  7)   872    251        1 ( 442/ 503- 572/ 679)  36.3 E50-E34 S19-S08
     14  275  601 (S12E64 L=343)   612     65        1 ( 259/ 567- 296/ 641)  59.0 E65-E61 S14-S09
     15  361  605 (S12E53 L=353)   421     55        0 ( 336/ 581- 395/ 628) -45.0 E55-E50 S13-S10
     16 1404  395 (S42W19 L= 66)   446     61        2 (1376/ 368-1443/ 415)   0.0 W16-W22 S44-S39
     17 1235  304 (S46W04 L= 50)   808    101        2 (1197/ 273-1293/ 334)  14.0 W00-W08 S48-S42
     18 1041  216 (S51E16 L= 30)  1154    136        2 ( 995/ 168-1087/ 276)  47.7 E22-E10 S56-S45
  

    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