Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1760 (N50E02 L=273)   432     71        3 (1314/1741-1378/1779)   9.5 E05-W01 N48-N53
      2  543 1676 (N20E52 L=224)   856    127        2 ( 500/1637- 601/1697)  45.0 E55-E47 N16-N21
      3  832 1482 (N14E23 L=252)   723    294        0 ( 790/1374- 874/1618) -47.3 E30-E17 N08-N21
      4  204 1235 (S13E62 L=213)  4476    731       21 ( 118/1104- 321/1371)  80.0 E72-E54 S20-S03
      5  692  998 (S17E19 L=257)   805    275        1 ( 573/ 961- 843/1028)  16.4 E26-E09 S19-S13
      6  210  943 (S30E55 L=220)  1577    278        2 ( 130/ 891- 272/ 986)  -7.6 E63-E48 S32-S27
      7  367  944 (S27E41 L=235)   450    127        0 ( 320/ 930- 420/ 961)  20.6 E44-E36 S28-S25
      8 1361  849 (S08W24 L=299)   717    191        1 (1308/ 759-1401/ 930) -35.0 W18-W28 S11-S04
      9  813  815 (S24E07 L=269)  1941    713        0 ( 498/ 719-1137/ 897)  11.1 E28-W14 S33-S16
     10 1345  604 (S21W32 L=307)   399    107        0 (1307/ 557-1381/ 642) -23.2 W27-W35 S22-S19
     11  994  511 (S37W14 L=290)  1069    256        2 ( 867/ 457-1112/ 545)  15.5 W05-W24 S41-S34
     12 1600  414 (S22W60 L=335)   400     60        1 (1577/ 386-1622/ 441) -76.0 W57-W61 S24-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