Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1440 1735 (N32W55 L=101)   807    215        2 (1383/1699-1508/1756)  -9.5 W48-W60 N30-N34
      2 1008 1243 (N19W05 L= 52)   763    225        1 ( 914/1218-1143/1274) -30.3 W00-W12 N14-N22
      3 1512 1212 (N02W33 L= 80)   475    186        2 (1473/1182-1550/1239)   9.5 W30-W36 N00-N04
      4 1885 1155 (S13W63 L=110)  2748    448        6 (1831/1085-1949/1221)  15.9 W54-W72 S17-S08
      5 1676  807 (S25W35 L= 81)  2201    506        0 (1548/ 674-1734/1055)  48.2 W23-W43 S32-S12
      6 1040  892 (S01E03 L= 44)   339    119        0 (1020/ 836-1063/ 957)  63.4 E03-W00 S05-N01
      7  802  891 (N05E16 L= 31)   225     60        0 ( 768/ 873- 846/ 909) -36.9 E17-E13 N02-N06
      8 1202  809 (S10W04 L= 51)   416    161        0 (1119/ 775-1274/ 858)   0.0 E01-W09 S11-S09
      9 1255  687 (S19W04 L= 51)   750    241        1 (1224/ 577-1315/ 809)  47.7 W00-W10 S25-S13
     10  290  672 (N03E59 L=348)  1346    342        1 ( 203/ 624- 375/ 745)   2.7 E69-E48 N00-N06
     11  883  478 (S21E24 L= 23)   256     63        2 ( 862/ 451- 902/ 508)  21.8 E26-E21 S23-S19
  

    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