Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1455 1701 (N45W30 L= 62)   738    142        1 (1408/1674-1507/1729)  20.0 W24-W35 N41-N47
      2 1534 1652 (N41W36 L= 68)   451     94        0 (1515/1619-1551/1684)  50.2 W32-W38 N38-N44
      3 1289 1593 (N33W14 L= 46)   693    123        0 (1237/1531-1327/1640) -54.5 W10-W16 N29-N36
      4 1645 1606 (N39W45 L= 77)   460     96        0 (1622/1580-1672/1627)  29.7 W41-W48 N36-N41
      5 1720 1576 (N38W52 L= 84)   583    113        1 (1699/1537-1745/1602)  45.0 W48-W55 N34-N40
      6  411 1368 (N13E44 L=349)   704    179        1 ( 345/1321- 467/1403) -15.3 E49-E38 N09-N14
      7 1246 1256 (N11W12 L= 44)   259     54        1 (1215/1232-1285/1280) -11.3 W09-W14 N09-N12
      8 1783 1210 (N13W51 L= 83)   381     74        0 (1760/1167-1800/1246) -45.0 W48-W52 N11-N15
      9 1560 1132 (N06W33 L= 65)   708    237        2 (1511/1001-1626/1244) -58.4 W29-W38 S01-N12
     10  185  768 (S24E62 L=330)  1537    276        2 ( 160/ 658- 215/ 868)  65.2 E64-E57 S30-S17
     11 1003  615 (S30W02 L= 34)   518    180        0 ( 924/ 556-1081/ 672) -31.0 E03-W07 S33-S26
     12  696  564 (S36E20 L= 12)   427     74        1 ( 647/ 538- 741/ 588)  -8.1 E23-E16 S37-S34
  

    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