Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1550 1712 (N30W66 L= 24)   470     61        1 (1516/1689-1579/1735)   0.0 W61-W72 N29-N31
      2 1201 1649 (N38W33 L=351)   388     74        1 (1176/1614-1230/1672)  76.0 W30-W35 N36-N40
      3 1590 1597 (N23W58 L= 16)   898    227        0 (1486/1565-1683/1618) -31.0 W50-W66 N18-N27
      4  621 1373 (N37E18 L=300)   329     75        1 ( 597/1336- 651/1399)  36.9 E20-E15 N34-N38
      5 1644 1306 (N06W46 L=  4)   625    134        1 (1617/1241-1668/1367)  66.8 W42-W48 N03-N10
      6  173  862 (N14E67 L=251)  2148    351        4 ( 118/ 726- 225/1028) -90.0 E75-E61 N05-N23
      7 1554  920 (S12W29 L=347)   711    185        1 (1526/ 844-1574/1016)  53.1 W25-W32 S15-S06
      8  914  755 (S05E13 L=306)   470    170        1 ( 892/ 677- 941/ 836)  60.9 E14-E09 S10-S01
      9  422  611 (S04E49 L=269)   312     52        0 ( 403/ 586- 448/ 634)  21.8 E51-E46 S05-S03
     10 1226  454 (S32E02 L=316)  1563    321        3 (1093/ 356-1318/ 556)  19.3 E14-W06 S36-S27
     11  719  459 (S19E34 L=284)   560     83        1 ( 696/ 417- 748/ 503)  63.4 E35-E31 S22-S16
     12  918  433 (S25E22 L=296)   411     97        0 ( 891/ 396- 948/ 476)  29.7 E25-E18 S27-S23
     13 1456  299 (S55W13 L=331)   535     83        2 (1417/ 279-1490/ 316)  -7.1 W08-W16 S57-S51
  

    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