Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1629 (N35E18 L=179)  3006    663        7 ( 718/1512-1214/1773)  -6.5 E44-W00 N30-N39
      2  391 1592 (N20E60 L=137)   963    316        1 ( 374/1504- 410/1679) -45.0 E67-E55 N14-N25
      3 1286 1498 (N33W07 L=204)   289     72        0 (1259/1475-1327/1514)  36.9 W04-W09 N30-N34
      4  613 1453 (N15E36 L=162)   383     74        0 ( 590/1418- 644/1496)  80.5 E36-E34 N12-N18
      5  378 1276 (N00E48 L=150)   516     84        2 ( 340/1247- 414/1311)  45.0 E50-E45 S01-N02
      6 1592 1078 (N13W35 L=232)   271     51        0 (1574/1044-1612/1115) -63.4 W32-W36 N11-N15
      7 1857 1004 (N15W60 L=257)   778    156        1 (1815/ 972-1889/1036) -11.3 W54-W64 N12-N16
      8  754  963 (S11E15 L=182)   229     61        0 ( 723/ 935- 786/ 997) -21.8 E17-E12 S11-S09
      9  918  912 (S10E04 L=193)   345     82        1 ( 879/ 887- 963/ 934)  31.0 E06-E01 S11-S08
     10 1667  684 (S08W51 L=248)   694    156        0 (1652/ 623-1693/ 750) -49.4 W47-W55 S11-S04
     11 1326  626 (S19W28 L=225)   264     65        0 (1291/ 612-1361/ 639)  36.9 W25-W30 S20-S17
     12 1011  571 (S29W09 L=207)   379     86        1 ( 982/ 550-1053/ 601)  21.8 W07-W12 S30-S27
  

    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