Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1395 1799 (N39W60 L=287)  2388    245        6 (1289/1763-1519/1830) -12.5 W46-W74 N34-N42
      2 1655 1614 (N20W65 L=292)  2623    312        3 (1559/1520-1782/1720) -67.6 W60-W73 N11-N28
      3 1281 1614 (N31W35 L=262)   676    205        1 (1188/1585-1379/1633) -32.5 W28-W39 N26-N33
      4  990 1447 (N29W09 L=236)   219     93        1 ( 948/1425-1032/1470) -51.3 W07-W11 N26-N31
      5  885 1435 (N31W02 L=229)   213     72        0 ( 846/1414- 916/1462)  14.0 W00-W04 N29-N31
      6  598 1413 (N37E19 L=209)   552    230        1 ( 509/1385- 706/1454)  -3.4 E26-E09 N35-N38
      7 1265 1286 (N13W20 L=247)   470    144        0 (1225/1221-1307/1346) -90.0 W19-W21 N08-N16
      8 1771 1107 (S10W50 L=277)   377     84        0 (1743/1069-1793/1147) -90.0 W48-W52 S12-S07
      9  202  898 (N14E61 L=166)   843    101        1 ( 160/ 867- 248/ 931) -29.1 E66-E55 N11-N16
     10  905  815 (S05E12 L=216)   464    172        2 ( 854/ 736- 946/ 869) -77.5 E13-E10 S10-S01
     11 1831  768 (S33W54 L=281)   790    117        2 (1809/ 727-1868/ 831) -45.0 W50-W61 S37-S27
     12 1686  705 (S32W37 L=264)   701    131        0 (1636/ 670-1730/ 744)   6.3 W32-W41 S34-S29
     13 1175  632 (S22E01 L=226)   880    191        1 (1122/ 557-1209/ 724)  45.0 E05-W03 S26-S17
     14 1574  514 (S42W26 L=253)  1305    228        1 (1478/ 452-1640/ 593)  17.4 W16-W32 S45-S38
     15  535  375 (S23E53 L=174)   516     97        0 ( 520/ 328- 555/ 428)  33.7 E57-E48 S26-S20
     16  664  325 (S30E47 L=181)   343     62        0 ( 646/ 293- 681/ 362)  26.6 E50-E42 S31-S27
     17 1141  246 (S50E18 L=209)  1706    273        6 (1015/ 192-1268/ 302)   0.0 E32-E04 S51-S47
  

    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