Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1391 (N27W07 L=120)   404    165        1 ( 891/1362-1054/1433)  -5.2 W02-W13 N25-N28
      2  721 1396 (N34E08 L=104)   692    175        2 ( 646/1373- 811/1416) -28.6 E14-E02 N30-N36
      3 1374 1377 (N15W30 L=143)   485    185        0 (1293/1362-1462/1397) -24.4 W25-W36 N12-N17
      4 1805 1318 (S02W61 L=173)  1252    192        4 (1759/1279-1853/1355) -35.5 W56-W66 S04-N01
      5 1190 1266 (N14W16 L=128)   439     87        1 (1153/1245-1243/1296) -45.0 W13-W18 N12-N16
      6 1773 1253 (S04W55 L=167)   359    111        0 (1746/1232-1809/1268) -15.9 W51-W58 S05-S02
      7 1734 1065 (S13W45 L=157)  1165    345        3 (1713/ 914-1775/1232)  71.6 W41-W50 S21-S03
      8  187  948 (N20E60 L= 53)   484     88        1 ( 163/ 921- 209/ 984) -53.1 E61-E56 N18-N22
      9 1517  793 (S21W23 L=135)   617    176        2 (1479/ 703-1548/ 865)  84.8 W21-W24 S27-S16
     10 1134  491 (S29E09 L=103)   383    114        0 (1119/ 439-1148/ 541)  60.3 E11-E07 S31-S24
     11  611  359 (S23E48 L= 65)   869    182        2 ( 587/ 288- 639/ 433)  55.0 E51-E42 S28-S18
     12 1339  338 (S47E01 L=111)  1296    191        2 (1291/ 274-1391/ 381) -90.0 E03-W04 S53-S41
     13 1258  287 (S48E09 L=103)  1343    162        2 (1189/ 242-1313/ 349)   0.0 E16-E03 S51-S42
  

    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