Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  776 1813 (N56E40 L=307)  1622    117        2 ( 722/1751- 823/1851) -45.0 E48-E32 N50-N60
      2  988 1593 (N39E08 L=339)   359     85        1 ( 959/1573-1038/1609)  -8.1 E10-E03 N37-N40
      3  543 1488 (N28E40 L=307)   824    198        0 ( 465/1449- 634/1534) -13.2 E48-E31 N25-N30
      4 1528 1399 (N29W35 L= 21)   828    247        0 (1466/1309-1598/1480) -57.5 W30-W39 N23-N34
      5  336 1382 (N19E55 L=291)  1515    186        3 ( 277/1296- 377/1457) -27.9 E66-E48 N14-N23
      6 1290 1431 (N30W16 L=  2)   220     50        0 (1261/1418-1319/1447)   0.0 W13-W18 N28-N30
      7  797 1321 (N19E17 L=329)   573    206        0 ( 714/1257- 867/1369) -26.6 E23-E11 N15-N21
      8 1364 1071 (N07W21 L=  7)   449    100        0 (1339/1032-1392/1118)  63.4 W19-W22 N04-N10
      9 1610  554 (S24W48 L= 34)   436    105        0 (1567/ 531-1658/ 590) -14.0 W42-W54 S25-S22
     10 1473  355 (S40W47 L= 33)   561     72        0 (1445/ 328-1492/ 395) -71.6 W43-W50 S42-S36
     11 1118  315 (S46W15 L=  1)   312     48        1 (1095/ 298-1143/ 338) -31.0 W12-W17 S47-S44
     12 1270  296 (S47W30 L= 16)   495     29        0 (1244/ 273-1294/ 318) -23.2 W26-W33 S48-S44
  

    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