Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1677 1565 (N26W61 L= 81)   468     83        0 (1644/1525-1706/1595) -68.2 W59-W62 N23-N28
      2  717 1342 (N30E16 L=  5)  1319    294        2 ( 578/1298- 854/1397)   2.7 E26-E05 N27-N31
      3 1345 1235 (N14W23 L= 44)   642    103        3 (1286/1198-1398/1269) -29.7 W19-W26 N11-N15
      4  201 1041 (N16E61 L=320)   890    230        0 ( 151/1005- 259/1085) -22.6 E66-E54 N13-N18
      5 1693  980 (S08W43 L= 64)   849    255        0 (1639/ 839-1735/1084)  49.8 W37-W48 S15-S02
      6  715  744 (S07E23 L=358)  1358    435        2 ( 641/ 591- 759/ 893)  59.5 E31-E19 S15-N02
      7  257  798 (N01E59 L=322)   522     76        0 ( 244/ 752- 274/ 837) -90.0 E60-E56 S02-N03
      8  980  809 (S06E06 L= 15)   458    129        2 ( 910/ 782-1043/ 831) -23.2 E10-E02 S08-S04
      9 1774  756 (S24W50 L= 71)  1117    160        3 (1719/ 723-1814/ 792)   6.3 W44-W53 S25-S21
     10 1599  710 (S23W34 L= 55)   375     85        0 (1586/ 670-1612/ 752) -71.6 W32-W34 S25-S19
     11  885  481 (S26E18 L=  3)   628    150        1 ( 829/ 424- 940/ 541)  35.0 E23-E13 S29-S22
     12  676  302 (S37E41 L=340)  1257    229        2 ( 637/ 253- 714/ 383)  57.5 E46-E37 S41-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