Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  533 1584 (N15E46 L=141)   334    126        0 ( 488/1567- 580/1604)   0.0 E50-E41 N14-N16
      2  293 1479 (N05E63 L=124)  1843    235        4 ( 224/1408- 354/1570)  65.6 E69-E57 N00-N11
      3 1662 1303 (N26W34 L=220)   470     90        1 (1623/1279-1698/1329)  33.7 W30-W36 N22-N27
      4 1319 1126 (N05W14 L=201)   304    107        0 (1258/1101-1355/1164)  56.3 W11-W15 N02-N08
      5  826 1050 (S10E12 L=175)   294    107        0 ( 809/ 995- 840/1102) -63.4 E13-E10 S13-S07
      6  607  925 (S22E23 L=164)   582    232        2 ( 492/ 862- 703/ 986)  37.6 E30-E17 S27-S16
      7 1561  481 (S23W51 L=238)   713    162        0 (1482/ 429-1607/ 542)  71.6 W48-W53 S27-S18
      8 1391  504 (S27W38 L=225)   479     94        2 (1336/ 477-1446/ 535)  50.2 W35-W40 S29-S23
      9  863  209 (S65W23 L=210)   555     62        1 ( 836/ 187- 885/ 238)   7.1 W19-W28 S67-S62
  

    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