Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  453 1667 (N59E55 L= 29)   686     71        0 ( 418/1643- 495/1685)   9.5 E63-E48 N56-N60
      2  902 1593 (N45W02 L= 86)   283     65        2 ( 865/1576- 939/1603) -11.3 E01-W04 N43-N46
      3 1600 1416 (N19W48 L=132)   655    230        0 (1545/1362-1664/1449) -45.0 W42-W51 N14-N21
      4 1322 1222 (N13W22 L=106)   232     81        1 (1280/1203-1358/1243)   0.0 W18-W24 N12-N14
      5  267 1062 (N17E55 L= 30)   571     62        0 ( 240/1030- 294/1086) -31.0 E57-E51 N15-N19
      6 1863 1003 (S11W62 L=147)   726     72        1 (1836/ 978-1897/1035)   0.0 W58-W67 S12-S08
      7 1038  569 (S22E06 L= 78)   951    242        0 ( 948/ 493-1120/ 661)  31.6 E13-W00 S26-S17
      8 1506  553 (S32W27 L=111)   400    129        0 (1465/ 499-1549/ 602)  45.0 W23-W30 S35-S29
      9  560  567 (S16E39 L= 46)   444     90        1 ( 530/ 525- 601/ 601)  20.6 E42-E34 S18-S13
     10 1093  464 (S30E05 L= 80)   879    301        0 ( 973/ 394-1213/ 531)  18.4 E14-W04 S33-S27
     11 1557  501 (S38W32 L=116)   446     76        1 (1519/ 482-1591/ 530)  18.4 W28-W35 S38-S35
     12  667  425 (S27E36 L= 48)   583     86        1 ( 645/ 387- 694/ 465)  71.6 E38-E33 S29-S23
     13 1343  442 (S37W13 L= 98)   284     67        0 (1311/ 418-1380/ 462)   9.5 W10-W16 S38-S35
  

    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