Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1755 1361 (N18W54 L=276)  1412    316        0 (1667/1240-1841/1478) -51.3 W48-W60 N10-N25
      2  328 1268 (N14E48 L=174)   758    187        2 ( 259/1228- 403/1290)  -8.1 E54-E40 N10-N15
      3 1638 1201 (N07W41 L=263)   243     51        0 (1619/1183-1660/1224) -33.7 W38-W42 N06-N08
      4  155 1097 (N04E65 L=157)  3347    381        6 (  94/ 976- 224/1206)  59.7 E76-E56 S02-N11
      5  483  865 (S11E35 L=187)  1375    422        2 ( 431/ 686- 524/1078)  87.5 E39-E33 S22-N01
      6 1042  864 (S13W01 L=223)   274    105        1 ( 983/ 837-1072/ 893) -36.9 E02-W02 S14-S10
      7 1840  846 (S13W60 L=283)   374     29        1 (1817/ 827-1857/ 863) -11.3 W57-W63 S14-S12
      8 1017  791 (S17E01 L=221)   268    131        0 ( 958/ 776-1077/ 821) -18.4 E04-W02 S17-S15
      9  868  655 (S25E11 L=211)   201    117        0 ( 829/ 623- 922/ 691)  33.7 E13-E07 S27-S23
     10 1697  564 (S32W54 L=276)   728    155        1 (1657/ 505-1752/ 626)  63.4 W51-W57 S35-S27
     11 1549  456 (S39W44 L=266)  1151    228        1 (1444/ 433-1640/ 482)   2.9 W32-W53 S41-S37
  

    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