Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  468 1512 (N18E44 L=  4)   788    121        2 ( 417/1469- 517/1565)  45.0 E47-E41 N14-N21
      2  311 1463 (N14E57 L=351)  4663    390       15 ( 224/1320- 455/1563)  53.1 E65-E46 N04-N20
      3 1716 1500 (N35W47 L= 95)   448     43        1 (1688/1479-1743/1525) -45.0 W43-W49 N33-N36
      4  525 1462 (N15E38 L= 10)   234     78        2 ( 496/1429- 546/1488) -76.0 E40-E35 N12-N16
      5  239 1034 (S12E54 L=354)  1558    279        5 ( 192/ 940- 291/1159)  77.0 E59-E49 S18-S05
      6  156 1084 (S10E64 L=344)   832     82        3 ( 133/1035- 174/1131)  80.5 E67-E61 S12-S06
      7 1848 1073 (N10W56 L=104)   801    126        0 (1826/1007-1867/1130)  82.9 W54-W58 N05-N13
      8  987  840 (S17W00 L= 48)   686    201        0 ( 947/ 749-1014/ 935) -70.0 E03-W02 S22-S11
      9 1532  766 (S14W36 L= 84)  1278    453        0 (1449/ 616-1610/ 929)  90.0 W31-W39 S22-S03
     10  284  851 (S23E51 L=357)  2925    546        3 (  96/ 774- 560/ 901)   8.1 E83-E27 S28-S18
     11 1212  853 (S13W13 L= 61)   259    115        0 (1162/ 838-1274/ 872)   0.0 W10-W17 S14-S12
     12  581  706 (S30E26 L= 22)  2263    461        6 ( 370/ 627- 728/ 851)  30.4 E45-E15 S36-S19
     13 1297  755 (S18W20 L= 68)   530    269        0 (1200/ 662-1379/ 845)  61.7 W16-W23 S24-S11
     14 1785  747 (S10W57 L=105)   306     70        0 (1749/ 731-1810/ 771)  31.0 W53-W59 S11-S08
     15 1125  406 (S43W18 L= 66)  2494    442        7 ( 920/ 368-1335/ 455)   8.6 W01-W35 S45-S38
     16  812  391 (S49E07 L= 41)   243     36        0 ( 792/ 373- 834/ 410)  45.0 E08-E05 S50-S47
     17 1596  369 (S36W62 L=110)   705     79        0 (1559/ 351-1634/ 392)  12.5 W56-W66 S37-S34
  

    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