Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  416 1266 (N17E43 L=220)   823    137        2 ( 373/1214- 447/1336)  90.0 E46-E39 N13-N21
      2  925 1200 (N14E07 L=256)   512    228        3 ( 873/1130-1012/1296) -45.0 E10-W00 N09-N20
      3  938 1280 (N19E06 L=257)   220     52        0 ( 910/1266- 969/1294)   0.0 E07-E03 N18-N20
      4 1097 1196 (N14W04 L=267)   528    221        1 (1029/1137-1151/1288)  60.9 W00-W07 N10-N19
      5  212 1192 (N11E62 L=200)   881    180        1 ( 169/1138- 241/1235) -40.6 E68-E58 N07-N13
      6  266 1147 (N09E55 L=207)   367     60        0 ( 248/1118- 283/1178)  76.0 E57-E53 N06-N10
      7 1244 1111 (N09W14 L=276)   255     68        1 (1215/1079-1268/1150)  68.2 W11-W15 N06-N11
      8  991  862 (S07E02 L=261)   253     97        0 ( 961/ 829-1039/ 880)  26.6 E03-W01 S08-S05
      9 1072  775 (S12W03 L=266)   291     28        1 (1046/ 746-1095/ 802) -33.7 W01-W04 S13-S10
     10  306  705 (S19E53 L=210)  1294    241        4 ( 246/ 627- 353/ 796) -65.6 E57-E48 S24-S13
     11 1396  743 (S14W25 L=287)   548     95        1 (1368/ 700-1426/ 790)  59.0 W22-W26 S16-S11
     12 1106  716 (S16W06 L=268)   200     79        0 (1078/ 682-1136/ 756)  51.3 W03-W07 S18-S13
     13 1000  608 (S23E01 L=261)   505    183        1 ( 980/ 516-1023/ 672)  84.8 E02-W00 S29-S18
     14  271  641 (S24E58 L=205)   656    124        1 ( 240/ 618- 302/ 669) -18.4 E61-E54 S25-S21
     15  504  596 (S25E37 L=226)   481    168        1 ( 453/ 557- 545/ 655)  33.7 E41-E32 S28-S21
     16  706  559 (S27E22 L=241)   421    127        0 ( 688/ 506- 729/ 607)  76.0 E23-E19 S31-S23
     17 1267  559 (S26W17 L=280)   633    136        1 (1195/ 535-1310/ 597)  14.0 W12-W20 S28-S23
     18 1058  487 (S32W03 L=266)   262     52        0 (1035/ 458-1078/ 516)  53.1 W01-W04 S33-S29
     19  402  431 (S39E53 L=210)   553    148        2 ( 386/ 403- 419/ 462)  45.0 E56-E49 S41-S36
  

    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