Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1573 (N29W61 L= 60)   888    183        0 (1599/1507-1691/1613)  10.0 W53-W70 N25-N30
      2 1526 1516 (N28W46 L= 45)   339     53        0 (1490/1501-1555/1530) -11.3 W43-W48 N26-N28
      3 1734 1454 (N19W63 L= 62)  2106    188        6 (1677/1396-1799/1530) -73.3 W59-W69 N15-N25
      4 1347 1459 (N28W30 L= 28)   345    123        0 (1329/1405-1383/1512) -90.0 W28-W31 N23-N31
      5  596 1336 (N31E26 L=333)  1406    296        6 ( 553/1215- 671/1460) -62.1 E30-E20 N22-N39
      6  863 1184 (N18E08 L=351)   238     49        2 ( 835/1164- 890/1205) -45.0 E09-E06 N16-N19
      7 1545 1107 (N03W34 L= 33)   502     74        1 (1511/1078-1573/1148)  36.9 W31-W36 N01-N05
      8 1721 1032 (S05W47 L= 46)   427    109        0 (1707/ 972-1735/1077) -90.0 W45-W48 S08-S01
      9 1594  930 (S08W35 L= 34)   239     60        0 (1567/ 906-1613/ 961)  36.9 W33-W37 S09-S06
     10  275  710 (S06E60 L=299)   691    104        1 ( 252/ 665- 297/ 761)  36.9 E63-E55 S08-S02
     11  971  632 (S17E09 L=350)   337    157        1 ( 910/ 594-1037/ 686)  16.7 E13-E03 S19-S14
  

    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