Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  774 1822 (N65W11 L= 40)   502     39        0 ( 757/1788- 796/1846)  36.9 W07-W15 N61-N68
      2 1101 1820 (N51W39 L= 67)   427     57        0 (1074/1788-1129/1845)  10.3 W33-W44 N48-N52
      3  685 1799 (N67E02 L= 27)  1240    105        1 ( 645/1747- 722/1841)  28.3 E08-W06 N61-N71
      4 1506 1703 (N30W59 L= 87)   759    105        1 (1474/1670-1551/1733)   5.2 W53-W65 N28-N30
      5  780 1712 (N54W04 L= 33)   408     54        0 ( 741/1696- 811/1730) -15.9 W00-W07 N51-N55
      6  938 1647 (N43W14 L= 42)   380     78        1 ( 915/1609- 962/1690)  45.0 W10-W15 N40-N46
      7  732 1651 (N50E02 L= 26)   404     59        1 ( 700/1628- 758/1676)   0.0 E05-W00 N48-N52
      8 1199 1527 (N28W25 L= 53)   466    189        0 (1148/1489-1280/1590) -60.9 W23-W29 N24-N33
      9 1259 1240 (N10W19 L= 47)   987    202        2 (1203/1151-1310/1330)  48.4 W13-W22 N06-N15
     10  604 1294 (N29E20 L=  8)   471    167        0 ( 529/1240- 661/1329) -45.0 E25-E16 N24-N33
     11  942 1272 (N20W01 L= 30)   341    126        0 ( 879/1250-1003/1291) -39.8 E02-W04 N17-N22
     12  145  945 (N18E68 L=320)  5094    385        5 ( 101/ 824- 200/1126) -83.7 E78-E61 N10-N28
     13  203 1081 (N25E58 L=330)   446     54        0 ( 184/1054- 220/1111)  45.0 E61-E56 N23-N26
     14 1756  857 (S25W45 L= 73)   486     82        0 (1709/ 845-1804/ 869) -18.4 W40-W49 S26-S22
     15 1008  306 (S40E24 L=  5)  1087    256        1 ( 926/ 224-1093/ 410)  22.2 E35-E13 S45-S34
     16 1470  297 (S59W15 L= 43)   510     43        0 (1447/ 279-1497/ 320) -50.2 W12-W18 S61-S55
     17 1309  212 (S61E07 L= 22)   636     63        1 (1272/ 196-1350/ 234) -18.4 E10-E01 S63-S57
  

    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