Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  627 1632 (N28E42 L=113)   471    110        2 ( 598/1600- 648/1681) -63.4 E44-E39 N25-N31
      2 1432 1486 (N36W19 L=174)  3501    813       15 (1111/1344-1720/1632)  -4.6 E07-W44 N30-N39
      3  794 1509 (N23E25 L=130)   355    154        0 ( 754/1452- 849/1557)   6.3 E29-E19 N20-N25
      4  617 1437 (N15E35 L=120)   346    155        3 ( 565/1403- 654/1462)   0.0 E38-E31 N13-N16
      5 1251 1265 (N17W09 L=164)   542    223        2 (1189/1205-1347/1299)  26.6 W05-W15 N12-N19
      6 1253 1121 (N09W12 L=167)   280    101        2 (1200/1102-1306/1145)   9.5 W08-W14 N07-N10
      7  970 1044 (S01E04 L=151)   329     96        0 ( 925/1016-1016/1082)  45.0 E05-E01 S03-N01
      8  372  753 (S30E41 L=114)   322     90        0 ( 351/ 722- 401/ 778)  38.7 E42-E37 S32-S28
      9  999  757 (S18W04 L=159)   224     34        1 ( 976/ 740-1027/ 776)   0.0 W02-W05 S18-S16
     10 1528  712 (S10W38 L=193)  1263    309        6 (1410/ 678-1629/ 766)  18.4 W29-W45 S12-S04
     11 1759  648 (S09W61 L=216)   745    128        1 (1724/ 592-1786/ 688)  90.0 W57-W63 S12-S05
     12 1416  408 (S31W43 L=198)   591    170        2 (1359/ 390-1477/ 441)   0.0 W37-W47 S32-S29
     13 1138  272 (S48W32 L=187)   335     47        0 (1114/ 254-1162/ 292)  76.0 W29-W33 S49-S45
  

    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