Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1397 1633 (N44W32 L=130)  1371    275        4 (1341/1567-1434/1693)  34.7 W25-W38 N38-N48
      2  927 1476 (N32E07 L= 91)   335     70        1 ( 910/1436- 950/1518) -80.5 E08-E05 N29-N35
      3 1723 1459 (N30W56 L=154)   463     85        1 (1700/1435-1747/1478) -45.0 W53-W58 N28-N31
      4 1553 1311 (N21W37 L=135)   454     82        0 (1513/1272-1595/1346) -45.0 W33-W39 N18-N23
      5  748 1248 (N17E18 L= 80)   308     98        3 ( 724/1215- 783/1298)  59.0 E19-E16 N15-N20
      6 1551 1225 (N15W35 L=133)   561    174        2 (1481/1202-1614/1259)   5.7 W30-W40 N14-N17
      7 1095 1162 (N12W04 L=102)   365    110        0 (1076/1115-1124/1220) -71.6 W03-W06 N09-N15
      8  800 1181 (N13E14 L= 84)   221     70        0 ( 770/1148- 820/1215) -53.1 E16-E12 N10-N15
      9 1557 1096 (N07W35 L=133)   389    120        0 (1501/1076-1609/1128) -20.6 W31-W39 N06-N09
     10  492  997 (N01E35 L= 63)   365    118        0 ( 442/ 966- 545/1033)  20.6 E38-E30 N00-N03
     11 1356  783 (S12W21 L=119)   205     64        0 (1340/ 765-1384/ 813) -45.0 W20-W23 S12-S09
     12  387  701 (S18E45 L= 53)   300     65        0 ( 375/ 665- 399/ 732)  68.2 E46-E43 S20-S15
     13 1574  624 (S23W39 L=137)   591     80        2 (1534/ 588-1611/ 660)   9.5 W36-W42 S24-S20
     14  708  626 (S22E21 L= 77)   654    107        0 ( 652/ 601- 768/ 656)  -7.1 E25-E16 S23-S20
     15 1104  400 (S38W07 L=105)   407     92        1 (1082/ 361-1124/ 445) -81.9 W04-W07 S41-S34
     16  986  350 (S43E03 L= 95)  1179    259        3 ( 916/ 318-1059/ 414)  10.3 E08-W03 S45-S37
     17  712  375 (S41E26 L= 72)   419     77        0 ( 683/ 355- 734/ 409)  38.7 E28-E23 S42-S38
  

    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