Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1465 1535 (N25W44 L= 62)   706    209        0 (1378/1492-1564/1574) -42.0 W39-W49 N20-N29
      2 1248 1420 (N25W24 L= 41)   627    205        2 (1191/1326-1287/1499)  33.7 W17-W29 N20-N28
      3  903 1320 (N27E01 L= 17)   984    326        0 ( 814/1192- 976/1441) -80.5 E04-W01 N18-N36
      4 1833 1305 (N02W66 L= 84)  1327    197        2 (1778/1220-1883/1374) -65.6 W62-W70 S04-N07
      5  783 1139 (N19E12 L=  5)   288    109        1 ( 721/1111- 823/1166) -33.7 E15-E09 N16-N21
      6  224  969 (N18E60 L=318)   945    241        0 ( 209/ 852- 258/1058)  69.0 E62-E54 N10-N23
      7 1785  957 (S16W50 L= 67)   415    100        0 (1739/ 939-1819/ 986)   0.0 W45-W54 S17-S14
      8 1671  764 (S25W37 L= 54)   492    107        1 (1656/ 715-1686/ 818) -90.0 W35-W38 S28-S21
      9  835  731 (S06E17 L=360)   291    124        1 ( 805/ 682- 884/ 770)  33.7 E20-E14 S08-S04
     10 1108  512 (S26E07 L= 10)   456    193        0 (1034/ 462-1188/ 563)  13.0 E13-W00 S27-S24
     11  809  502 (S20E26 L=352)   614    232        0 ( 702/ 478- 876/ 556)   4.4 E33-E20 S21-S17
  

    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