Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1655 (N31E18 L=278)   530    210        0 ( 925/1585-1039/1729)  81.3 E20-E16 N24-N37
      2 1114 1637 (N33E08 L=288)   602    167        0 (1064/1568-1157/1703) -38.7 E13-E03 N28-N36
      3  525 1565 (N16E45 L=251)   993    180        1 ( 495/1517- 576/1653) -42.0 E51-E39 N13-N22
      4  311 1584 (N15E69 L=227)  3446    168        7 ( 253/1526- 364/1652)  60.9 E78-E59 N11-N20
      5 1535 1592 (N41W24 L=321)   854    155        1 (1457/1560-1601/1613)   8.7 W17-W30 N38-N43
      6 1769 1473 (N40W49 L=346)   551     92        0 (1717/1450-1807/1493)   4.4 W42-W55 N39-N41
      7 1165 1294 (N11W03 L=299)   628    197        1 (1130/1211-1211/1381)  85.2 W01-W04 N05-N17
      8 1329 1265 (N13W13 L=309)   385    127        0 (1310/1211-1345/1321) -63.4 W11-W14 N10-N16
      9 1524 1119 (N09W28 L=324)   543    182        0 (1474/1052-1562/1194)  73.3 W25-W29 N04-N14
     10 1203 1090 (N00W09 L=305)   436    167        1 (1127/1061-1275/1139)  -5.7 W03-W13 N00-N01
     11 1920  906 (N09W70 L=  6)  1419    111        3 (1886/ 859-1963/ 943)  50.2 W65-W81 N05-N11
     12 1028  907 (S14W03 L=299)   253    154        0 ( 979/ 860-1095/ 932)  45.0 W00-W06 S17-S11
     13  671  892 (S21E19 L=277)   246     86        0 ( 627/ 867- 707/ 918)  45.0 E21-E16 S23-S19
     14 1224  726 (S20W19 L=315)   450    144        2 (1171/ 676-1275/ 782)  63.4 W16-W20 S23-S15
     15 1464  654 (S18W36 L=333)   290    102        0 (1424/ 636-1502/ 688)  45.0 W33-W38 S19-S15
     16 1246  409 (S38W33 L=330)   362     80        0 (1211/ 393-1288/ 439)   9.5 W29-W37 S39-S36
     17 1569  387 (S30W61 L=357)   581    129        0 (1505/ 357-1616/ 421)  45.0 W56-W64 S32-S26
  

    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