Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1508 1591 (N41W19 L= 21)  1903    200        3 (1373/1476-1622/1680) -21.8 W07-W29 N35-N45
      2 1127 1642 (N33E10 L=352)   296     57        1 (1102/1616-1155/1670) -36.9 E11-E06 N31-N34
      3  740 1520 (N16E29 L=333)   423    150        0 ( 717/1451- 759/1579) -53.1 E31-E25 N11-N19
      4 1179 1325 (N13W02 L=  4)   210     70        0 (1155/1296-1196/1363)  90.0 W00-W02 N11-N15
      5 1349 1275 (N15W13 L= 15)   380    152        0 (1329/1193-1375/1349) -60.9 W11-W16 N10-N19
      6  214 1250 (S09E62 L=300)   576     91        4 ( 189/1214- 235/1294) -63.4 E64-E59 S10-S06
      7  365 1114 (S15E44 L=318)   264     58        1 ( 346/1086- 391/1135)  26.6 E45-E41 S16-S13
      8  388 1019 (S20E41 L=321)   362    115        0 ( 368/ 952- 402/1073)  90.0 E42-E39 S24-S16
      9  902  975 (S13E06 L=356)   774    200        0 ( 812/ 882- 968/1065)  72.9 E09-E03 S19-S06
     10 1603 1036 (N08W34 L= 36)   257     60        1 (1572/1017-1634/1052)  26.6 W32-W36 N06-N08
     11 1513  818 (S07W34 L= 36)  1288    335        2 (1374/ 677-1632/ 972)  69.1 W29-W37 S17-N04
     12 1180  672 (S23W18 L= 20)   801    199        2 (1106/ 642-1277/ 736)   0.0 W12-W24 S25-S20
  

    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