Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1339 1805 (N43W60 L=226)   512     59        2 (1307/1777-1362/1832)   4.8 W54-W67 N41-N43
      2 1623 1622 (N22W65 L=231)  2197    367        2 (1523/1524-1784/1685) -43.4 W56-W77 N10-N27
      3 1386 1462 (N22W35 L=201)  1232    460        4 (1215/1372-1560/1509) -41.6 W25-W43 N11-N28
      4  888 1355 (N29W01 L=166)   551    251        2 ( 785/1288- 973/1432) -61.7 E04-W04 N23-N36
      5  339 1218 (N34E45 L=121)   791    214        2 ( 288/1177- 411/1255) -21.0 E50-E37 N30-N36
      6  325 1061 (N25E47 L=119)   879    239        0 ( 285/ 981- 388/1160) -63.4 E50-E41 N19-N31
      7  453 1087 (N24E35 L=131)   341    127        0 ( 414/1047- 492/1140) -54.5 E37-E32 N20-N27
      8  270  790 (N09E58 L=108)  1426    260        6 ( 222/ 680- 304/ 901)  81.9 E61-E53 N02-N16
      9  489  670 (S01E42 L=124)   493    103        0 ( 460/ 642- 529/ 712) -76.0 E43-E39 S03-N01
     10 1792  597 (S44W52 L=217)   594     50        1 (1771/ 579-1818/ 626) -32.5 W47-W58 S47-S40
     11  395  488 (S11E59 L=107)   483     65        1 ( 374/ 456- 423/ 526)  29.7 E63-E54 S13-S08
     12 1006  354 (S33E22 L=144)   290     51        0 ( 980/ 330-1024/ 380) -90.0 E22-E19 S35-S30
  

    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