Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  864 1748 (N53E19 L=158)   482     63        1 ( 831/1724- 897/1769)  31.0 E21-E15 N50-N54
      2  352 1538 (N33E61 L=115)  1933    194        3 ( 297/1486- 427/1613) -11.6 E84-E49 N29-N37
      3 1342 1552 (N38W24 L=200)   428    109        0 (1303/1532-1386/1571)  -8.1 W20-W27 N36-N39
      4 1415 1513 (N35W29 L=205)   383     47        0 (1389/1481-1439/1538)  38.7 W26-W31 N32-N36
      5  632 1458 (N29E30 L=147)  1127    231        2 ( 541/1424- 745/1491)  -6.7 E37-E20 N26-N31
      6 1448 1459 (N31W31 L=207)   276     49        1 (1424/1442-1468/1485)  14.0 W28-W32 N29-N32
      7 1461 1329 (N22W29 L=206)   368    101        0 (1444/1284-1473/1381)  71.6 W28-W30 N19-N25
      8  449 1137 (N08E38 L=138)   548    172        0 ( 407/1062- 484/1208) -52.1 E42-E35 N03-N12
      9 1790 1159 (N11W56 L=232)   598     96        0 (1747/1117-1821/1190) -35.5 W51-W58 N08-N13
     10 1889 1103 (N08W68 L=244)   761     99        1 (1869/1056-1918/1154) -56.3 W65-W72 N04-N10
     11  798 1047 (N04E14 L=162)   258     63        0 ( 778/1012- 826/1077)  56.3 E15-E12 N01-N05
     12 1843  972 (S00W61 L=237)   586     77        0 (1825/ 935-1859/1011)  76.0 W58-W62 S02-N02
     13 1760  752 (S14W54 L=231)  1067    123        2 (1717/ 704-1808/ 796) -21.8 W49-W60 S16-S10
     14  390  668 (S21E45 L=131)   360    131        0 ( 344/ 645- 448/ 715)  21.8 E50-E39 S22-S17
     15 1068  566 (S26W04 L=180)  1106    239        3 (1001/ 522-1156/ 649) -38.7 W00-W10 S28-S20
     16  876  570 (S26E09 L=167)   398    106        0 ( 825/ 546- 922/ 607)  33.7 E12-E05 S27-S23
     17 1370  214 (S56W45 L=221)   683     77        0 (1343/ 186-1389/ 253) -45.0 W39-W51 S59-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