Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1178 1745 (N53W01 L=218)   325     70        0 (1149/1729-1212/1756)  36.9 E02-W03 N50-N54
      2  446 1425 (N18E46 L=171)   586    209        0 ( 435/1338- 460/1521) -65.6 E49-E43 N12-N23
      3 1182 1399 (N26W05 L=222)   266    140        0 (1166/1342-1199/1468) -83.7 W04-W06 N21-N30
      4 1451 1414 (N31W25 L=242)   403     82        2 (1415/1370-1488/1442) -38.7 W21-W27 N27-N32
      5  272 1165 (S01E55 L=162)   376     70        0 ( 252/1136- 293/1195) -11.3 E57-E52 S02-N01
      6 1770  895 (N02W54 L=271)  5704    969       14 (1676/ 561-1849/1169) -74.9 W42-W66 S18-N19
      7  260  950 (S14E53 L=164)   440     96        0 ( 237/ 913- 286/ 981) -73.8 E55-E50 S15-S12
      8  212  835 (S22E58 L=159)  2375    391        5 ( 175/ 769- 292/ 914) -90.0 E62-E49 S26-S16
      9  701  758 (S20E17 L=200)   271    116        0 ( 651/ 722- 760/ 803) -20.6 E21-E13 S21-S17
     10 1497  692 (S14W36 L=253)  1512    428        5 (1362/ 538-1627/ 794)  64.7 W28-W44 S25-S06
     11 1546  650 (S16W40 L=258)   598    174        3 (1491/ 612-1608/ 683)  35.5 W37-W45 S18-S13
     12 1533  416 (S32W49 L=266)   369     61        1 (1507/ 400-1560/ 439)  14.0 W47-W51 S32-S29
  

    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